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Abstract 



We discuss radiative backreaction for global strings described by the Kalb-Ramond 
action with an analogous derivation to that for the point electron in classical electrody- 
namics. We show how local corrections to the equations of motion allow one to separate 
the self-field of the string from that of the radiation field. Modifications to this 'local 
backreaction approximation' circumvent the runaway solutions familiar from electrody- 
namics, allowing these corrections to be used to evolve string trajectories numerically. 
Comparisons are made with analytic and numerical radiation calculations from previous 
work and the merits and limitations of this approach are discussed. Finally, we comment 
on the relevance of this work to simulations of a cosmological network of global or axion 
strings. These methods can also be applied to describe gravitational, electromagnetic 
and other forms of string radiative backreaction. 

1 Introduction 

Global strings arise naturally in theories with a spontaneously broken global U{1) sym- 
metries, such as axion models, superstrings and some GUT models. The Goldstone 
bosons or axions, associated with the broken symmetry, can be radiated by an oscillat- 
ing global string. This radiation power can be calculated using perturbed trajectories, 
but the effect of radiative damping on the actual string motion remains to be adequately 
understood. This deficiency is a major obstacle impeding the quantitative description of 
string network evolution. Radiative backreaction is believed to determine the nature of 
long string small-scale structure and small loop creation sizes and, hence, the amplitude 
and spectrum of the resulting radiation background. 

It is appropriate to begin discussion of the problem of radiative backreaction by 
considering the point electron in classical relativistic electrodynamics [1,2,3]. In this case, 
the self and radiation fields of the electron can be distinguished easily since the self-field 
is of order whereas the radiation field fall-off is of order 1/R. Careful analysis of 

the equations of motion leads to the renormalisation of the electron mass Mj-en = f-^e 
by the Coulomb self-field and the first-order approximation to the radiation force — the 
Abraham-Lorentz force — is given by 
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where X^{t) is the position on the electron's worldhne at time r. However, the depen- 
dence of this force on leads to problems in numerical applications because (1) has 
exponentially increasing solutions. These unphysical 'runaway' solutions can only be 
suppressed by rewriting the equations of motion as an integro-differential equation. In 
the following we shall derive the analogue of (1) for a global string using an antisymmetric 
tensor formalism. 

2 Self-field renormalisation 

The essential features of global strings in flat space are exhibited in the simple U(l) 
Goldstone model, with action given by 

where $ is a complex scalar field which can be split into a massive (real) component 
(j) and a massless (real periodic) Goldstone boson d. The analytic treatment of global 
string dynamics is hampered by the topological coupling of the self field of the string 
to the Goldstone boson radiation field. However, we can exploit the well-known duality 
between a massless scalar field and a two-index antisymmetric tensor to replace 
the Goldstone boson ■& in (2) via the relation (fP'dfj,{} — \fs.^fxv\pd^B^P . Performing 
this transformation carefully and integrating over the massive degrees of freedom about 
the two-dimensional string worldsheet X'*(cr, r) [4], yields the flat-space Kalb-Ramond 
action [5,6,7], 

S=-iio j V^dadT -^J ^d^xH'' - 27T fa J B^/sV^dadr, (3) 

where i?^a/3 = dnBaf3+dpBfj_o,+daBp^ is the field strength of B^p, the worldsheet metric 
is 7a6 = g^lu^aX^'^bX'' with 7 = det(7a6), and Vap = e^-^daX^dbXp is the antisymmetric 
vertex operator. 

Varying the action (3) with respect to the worldsheet coordinates and the antisym- 
metric tensor yields the string equations of motion and the tensor field equations, 

r (4) 

^ 4^ ja/3 ^ ^TT U / dadr (5^ {x - X((7, t)) F'^^ . 

In the Lorentz gauge, d^B^'^ = 0, the field equations can be solved using standard 
Green's function methods in retarded time, 

Bap{x) = 27r/a J dadfDr^etix - X{a, f))Vaf3{^,f) , (5) 

where D^etix) = {1 /27r)9 {xo)d (x"^) is the retarded Green's function. Defining = x^ — 
Xfxip^f) and treating a — a{f) one can perform various standard manipulations [2,8,9] 
to obtain the Lienard-Wiechart potential and its derivative 



r=TB 

(6) 
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where 



-TR 



and tr < t. One can separate the radiation field from the self field 



by noting that the radiation field is a solution of the homogeneous wave equation. If one 
defines the advanced Green's function, DadvC^^) = il/27r)9{—xo)5{x'^), then by analogy 
with electrodynamics [2,8] one can define Green's functions for the self and radiation 
fields 2^6^(2;) = Dret(a^) + Di,dv{x) and 2D,-i,d{x) = D^et{x) - L'adv(a^)- Thus we can 
calculate the derivatives of the Lienard-Wiechart potentials for both components 




+ 



-TR 



A.X df V \A.X\ 
A.X df\ \A.X\ 



(7) 



where A'^lf^^.^y — and tr < t,T'j^> t. 

The equations of motion (4) are problematic because the field diverges as any point 
on the string is approached — > X^(cr, r), as for the point electron in classical electro- 
dynamics. The electron has a natural scale, the classical electron radius, which allows the 
renormalisation of the electron mass and the calculation of the first-order approximation 
to the radiation force (1). In a similar manner for strings, the short distance divergences 
are cut off by the natural scale of the string width 6. However, the string self-field also 
diverges at large distances and so we must introduce another arbitrary cut-off scale A. 
Formally, this scale corresponds to the integrals in (5), (6) and (7) running over a region 
of size A around the point in question, that is, we only consider the effects from neigh- 
bouring string points within a distance of A/2. Note that we define A in terms of an 
invariant length of string. Physically, it will correspond to a distance beyond which the 
long-range string fields become uncorrelated and begin to cancel, so we are essentially 
approximating the effect of the nearby string with a square or 'top hat' window function. 
Our expectation, then, is that an appropriate choice for A would be near the average 
curvature radius of the string. 

For fiat-space string dynamics, the conformal string gauge is usually employed. 
However, when considering problems in which the string energy decays it is best to use 
the temporal transverse gauge in which X^ — t = r and X.X' = with X^^ = (t, X). 
The equations of motion for the string (4) and the energy of the string become 



/io^X- -^^^^ ^ = f , //oe = /°, E = iJ,oJdae, 



(8) 



where = XV(1 - X'^) and = {f, f). 

We will now renormalize the equations of motion (8), though the details of the 
procedure will be described elsewhere [10] (see also refs. [8,9,11]). For the purposes of 
this paper, it suffices to note that the expression for the derivatives of the Lienard- 
Wiechart potentials of the self and radiation fields (7) can be expanded locally about 
some point on the string. As discussed above, this procedure strongly suggests a natural 
scale for the otherwise arbitrary renormalisation cut-off A near the average curvature 
radius of the string. After a detailed set of manipulations which involve use of the gauge 
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conditions, one can deduce that 



H 



self 



fa 



rj-rad 



2X4 

/a 



log(A/5) + 0(A2) , 
3e 



4--- 1 • • 



(9) 



' 2X X e 



A + 0(A2) 



where ^[;:iQ:/3] = ^^a/? + + Aapij.- Note that the self-field has no order A term. 

Ignoring terms of order A^, we can then deduce expressions for the self-force and the 
radiation backreaction force, 



f-'f = -27r/flog(A/5)(^X- 



f-d=7r/^A<'^eX + 



2e 



XX 



+ 3e 



+ 



[--( 


^X' 




3e ^ 







1 -X2, 

4(X-X)(X'-X) 
e (1-X2)2 



0,self 



-27r/flog(A/5)e. 
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(10) 

The expressions for f^^^^ and /O'^eif^ facilitate the well-known renormalisation of the 
equations of motion (8) and string energy. 



MA)(X-i[^) )=f-^ MA)6 = /°'-d 



E 



dae , 



(11) 



with the renormalised string tension /i(A) = fj,Q + 27r/f log(A/(5), while the expressions 
for f'^'^ and jO,rad i-gp^-esent the finite radiation backreaction force in what we denote as 
the 'local backreaction approximation'. The assumptions underlying (11) are that the 
dominant contribution to the integrals (5), (6) and (7) come from string segments close 
to the point under consideration and that the long-range fields are uncorrelated beyond 
the string radius of curvature. We believe these are reasonable assumptions for general 
physical situations, notably for strings in a realistic interacting network. 



3 Modified string equations of motion 

The equations of motion for the radiation force (10) and (11) are rather complicated. 
However, at least for the string solutions considered in ref. [11], that is, closed loops and 
long periodic strings, we can demonstrate that some of these terms are sub-dominant 
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(this has also been verified numericaUy) . In particular, it is possible to approximate 

r ■- 1 /X' ■ X \ .1 

(12) 



Tj-rad 



3X4 

.ad ^ 47r/f A 



eX 



e V 1 - X2 



X' 



/ 



0,rad ^ 47r/f A 



e^X-X 
1 -X2 



This simplified form of the equations of motion using (12) still has serious shortcom- 
ings because of the presence of the X term. The equations have unphysical, exponentially 
growing or 'runaway' solutions which will, for example, plague any potential numerical 
applications. Furthermore, one would be required to store information at three differ- 
ent timesteps, fundamentally changing the nature of a numerical algorithm. It appears, 
however, that both these problems can be circumvented by resubstituting the equations 
of motion, that is, making the approximations X = X"/e^ and X = X"/e^ in (12). The 
equations of motion then acquire an analogue of a viscosity term for which there are only 
damped solutions. 

The reason for the suppression of the exponentially growing solution becomes ap- 
parent if we consider simplified one-dimensional equations. 



X — X" = aX , 



X -X" ^ aX" , 



(13) 



where we have performed the resubstitution assuming that a is small. We now take an 
approximately periodic solution, X" —Q^X, and we substitute the ansatz X ~ e"^*. 
The solutions for (13) are given respectively by the roots of the following polynomials in 

m. 



/(m) = am'^ — rr? — fl^ , — * giiTi) — —m^ — aVt^m — Vt^ . (14) 

If we rewrite /(m) = {m^ + Am^Vt^ +B){am-C), then we see that A = aO^ + C(q!^), 
B = 0{a'^) and C — l + 0{a'^). If we ignore terms C(q;^), then the solutions of g{m) — 
are approximately solutions of /(m) = 0. However, the real positive solution of /(m) = 0, 
corresponding to the exponentially growing solution, is not a solution of g{m) = 0. 

It remains to recast the simplified resubstituted equations of motion into a first- 
order form accessible to numerical solution. Defining a = X' — eX and (3 = X.' + eX, 
the equations of motion can be rewritten as 



/.(A) 
/^(A)e 



a+ { - 



a 
e 



1 ( rad ^ f> 



O.rad 



X 



gjrrad _j_ yO,rad-j^ 



(15) 



/ 



o.rad 



X 



Using the above, one can then evolve string trajectories simply by modifying a total 
variation non-increasing (TVNI) algorithm [12,13] which has already been well-tested 
for string network evolution in an expanding universe. 
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Figure 1: Decay of S using the radiative backreaction force (dotted line) and numerical field 
theory simulations (solid line) for (a) a sinusoidal perturbation, (b) a helicoidal perturbation with 
unequal left- and right-moving amplitudes, and (c) a pure left moving helicoidal perturbation. Note 
the excellent quantitative agreement for all three cases. 

4 Numerical and analytic comparisons 

The power due to the radiation backreaction force can be found by differentiating the 
expression for the energy in (11) 



For closed loops of length L, we take the integration over the range < cr < L. Since 



recovers the well known result that the power loss from a loop is independent of its size 
L. 

More caution has to be exercised in comparisons with long string configurations 
because of the possibility of periodicity or other global correlations interfering and sup- 
pressing radiation power. The generic result for a long string solution parameterized by 
its wavelength L and amplitude to wavelength ratio e = 271 A/ L, where A is the ampli- 
tude, is X ~ 0{e), X ~ 0{eL~'^). Therefore the power per unit length is dP/dl — /3e'^/L. 
The efi^ect of this power loss can be calculated by analogy to the simple backreaction 




(16) 
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model presented in ref. [11]; it will lead to an exponential decay of the amplitude and 
oscillation energy per unit length, 



where a, (3 are constants dependent on the string configuration, defined as in ref. [11]. 

This result is in agreement with the exponential decay already shown to occur for 
general configurations with unequal left- and right- moving modes [11]. However, the 
periodic solutions discussed in ref. [14,11] were shown to have a weaker power law fall- 
off, dP/dl = Pe^/L, and some pure left-moving (or right-moving) solutions are known to 
propagate indefinitely along a straight string solution with no decay [15]. The resolution 
of this apparent contradiction lies in noticing that the renormalization procedure ignores 
any global correlations of the long-range fields outside the cut-off scale A. For example, 
the non-radiating left-moving solution has an accompanying left-moving perturbation 
in the long-range Goldstone field which is (unphysically) correlated out to infinity [15]. 
Accordingly, if we suppress the Goldstone field correlations at the curvature radius of 
the string perturbations, then we should recover the exponential decay anticipated in 



We have extensively tested our 'local backreaction approximation', using the modi- 
fied Nambu equations of motion (15), by comparing with field theory simulations of radi- 
ating strings in the Goldstone model (2) (see ref. [11]). Here, we shall only briefiy report 
on the remarkable correspondence of the radiative decay results, leaving further details 
for a longer publication [10]. Note that in the field theory simulations we evolve initial 
trajectories for which Goldstone field correlations are suppressed by gaussian smoothing 
for transverse distances greater than the string curvature radius. The massless field of a 
perturbed string approaches that of a straight string at large distances from the core, as 
we would expect in a general physical context for random string small-scale structure. 

Fig. 1 illustrates the excellent quantitative agreement for three different long string 
configurations, including a standing sinusoidal perturbation, a helicoidal perturbation 
with unequal left- and right- moving amplitudes, and a pure left-moving helicoidal per- 
turbation. In all three cases, the decay is seen to be exponential. The agreement persists 
for the longest time for the (generic) unequal left- and right-moving configuration be- 
cause exponential decay is predicted even after field correlations have relaxed at large 
distances. By comparison to the simple backreaction model for exponential decay (17) 
one can use the numerical field theory results to estimate A fa 0.1 L. We had anticipated 
that A should be normalized to a distance near the string radius of curvature R, which 
for a sinusoidal perturbation is i? ~ L/4. The fact that A is smaller than R validates 
taking only the lowest order terms in the expansion (10) for the radiative backreaction 
force. 

Fig. 2 compares the evolution of a sharp kink in the local backreaction approxima- 
tion with a field theory simulation. The results are almost indistinguishable except for 
the computational advantages of the former which, in this case, saved a factor of 10^ 
in cpu time and 10^ in allocated memory. Backreaction leads to a substantial rounding 
of the kink, in agreement with the intuitive picture described in ref. [11]. A more thor- 
ough investigation of this kink rounding, including spectral analysis, will be presented 




(17) 



(17). 



in ref. [10]. 
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Figure 2: Decay of kink perturbation (£o=0.9) for using (a) the radiation backreaction force 
and (b) numerical field theory. Notice the visible rounding of the kink in both cases 



5 Discussion and conclusions 

The local backreaction approximation which we have formulated appears to provide a 
good description of radiative damping in a realistic context, notably that of an evolving 
string network with uncorrelated left- and right- movers. In physical situations when the 
long-range string fields are only correlated out to distances comparable with the string 
curvature radius, numerical simulations confirmed the expected rate of exponential decay. 
This approach also appears to satisfactorily describe loop decay. Of course, our local 
approximation does not apply for solutions with global correlations in the long-range 
fields, but our expectation is that such solutions are not of physical relevance. These 
methods can also be applied to gravitational and electromagnetic backreaction which we 
will discuss, along with more detailed tests, in a longer publication on this subject [10]. 

These methods should prove useful for phenomeno logically describing radiative ef- 
fects in simulations of string networks. The inclusion of such backreaction terms may 
make tractable a number of outstanding cosmic string problems, enabling us to under- 
stand the true nature and existence of the 'scaling' solution, as well as the amplitude and 
spectrum of axion and gravitational wave backgrounds produced by radiating strings. 
Reformulated, this basic approach to describing string radiation might also find applica- 
tion in non-relativistic physical contexts, such as describing the radiation sound waves 
by vortex- lines in condensed matter systems. 
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